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We characterize a pair of Cooper-pair boxes coupled with a fixed capacitor using spectroscopy 
and measurements of the ground-state quantum capacitance. We use the extracted parameters to 
estimate the concurrence, or degree of entanglement between the two qubits. We also present a 
thorough demonstration of a multiplexed quantum capacitance measurement technique, which is in 
principle scalable to a large array of superconducting qubits. 



I. INTRODUCTION 

Entanglement is a unique feature of quantum mechan- 
ics which is of great fundamental interest, and is also an 
essential resource for quantum information processing. 
In the field of superconducting electronics, macroscopic 
entangled states of tw o or m or e qub its have r ecent ly been 
investigated in charge,'^^^flux,Sl3^and phas#'^ qubits . 
The pioneering experiments with fixed coupling explored 
the two-qubit state space through coherent oscillations^ 
and microwave spectroscopy,'^ mapping out the avoided 
level crossings characteristic of coupled qubits. Other 
key developments from the perspective of quantum infor- 
mation processing were the demonstration of controlled 
gate operationPI and simultaneous single-shot two-qubit 
measurement^^ in the case of fixed coupling. More 
recently, the strong-coupling limit of circuit QED has 
been used to exchange quantum informat ion between two 
qubits via a superconducting resonator,!^ and a vari- 
ety of tunable couplin g schemes have been demonstrated 
experiment ally.* ^ \^o \ in Purthermore, extensive work has 
been performed to investigate the ground state entangle- 
ment of systems of two, three, and four coupled qubits 
in the context of adiabatic quantum computing ,121121121 
and such measurements have been directly compared 
with spectroscopic characterization and Landau-Zener 
interferometry. 

In this work, we characterize a system composed of 
two single Cooper-pair box (SCB) charge qubits, coupled 
with a fixed capacitor, using a frequency-multiplexed 
quantum capacitance measurement technique. The qubit 
parameters, including the coupling energy between the 
two qubits, are extracted spectroscopically and via mea- 
surements of the ground state capacitance. These pa- 
rameters can then be used to estimate the degree of en- 
tanglement for the two-qubit system as a function of gate 
charge, although such an inference does not in itself con- 
stitute a direct observation of the entangled state, as can 
be established using state tomography.fSl 

We also present the multiplexed quantum capacitance 
measurement (QCM), which is an effective method for 
probing the state of multiple qubits with a single RF 



line. In a single qubit, the QCM is a dispersive measure- 
ment of the reactive response of an L C oscillator coupled 
capacitively to the qubit island.'^^'^ The capacitance of 
the qubit has, in addition to its geometric capacitance, a 
term which is determined by the second derivative of the 
qubit energy with respect to the gate charge, so that the 
overall capacitance of the oscillator depends on the qubit 
state. The state of the SCB can be obtained in a capac- 
itance measurement by monitoring the center frequency 
of the oscillator with RF reflectometry. The oscillator is 
tuned to a frequency much lower than the qubit energy 
level splitting, minimizing measurement back action and 
disturbance to the qubit. This also has the benefit of fil- 
tering high-frequency noise from the RF line. A primary 
advantage of the QCM technique as compared to some 
earlier measurement devices, such as the single-electron 
transistor, is its applicability for qubit state discrimina- 
tion at the degeneracy point, where dephasing due to low- 
frequency voltage fiuctuations is minimized. In princi- 
ple, the QCM technique can be used to perform quantum- 
limited measurements. 

As superconducting circuits grow in complexity, effi- 
cient multiplexing schemes are required for practical op- 
eration at millike Ivin temperatures. By coupling each 
qubit to a high-Q oscillator of a different frequency, an 
array of qubits can be efficiently read out with multi- 
plexed QCM by applying a frequency comb to a single 
RF line. In this work, we demonstrate the use of such a 
technique to characterize a system of two coupled qubits. 
A similar concept is also being used to read out arrays of 
superconducting radiation detectors .1^21 

In the four-level approximation, the Hamiltonian for a 
two-qub it sy stem coupled by a fixed capacitance Cm is 
given bjl^lSfil 
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where Ecj = 2(CsiCs2-'c'^ ) ^^"^ ^•^J ^'^^ charging and 
Josephson energies, respectively, for the j'th qubit, £".,„ = 
CsiCs2-c^ mutual coupling energy between the 

two qubits, and C^j is the total capacitance for the jth 
qubit. In Eqn. ([Tj), cr^ ^ are the standard Pauli matrices, 
Czi = (8) 1, crz2 = 1 fX" and CTzz = CTz (g) CTz. Note 
that the interqubit coupling is along the zz axis, and is 
not controllable with the gate voltages. The coupling 
between qubits introduces avoided level crossings and a 
gate voltage asymmetry into the energy spectrum. 

In section II, we discuss the multiplexed QCM mea- 
surement, while in the remainder of the paper we present 
a characterization of the entangled two-qubit system. 
Section Ila provides a short review of the QCM concept 
for a single qubit, while section lib presents the details 
of the experiment and a discussion of device fabrication. 
In section lie, we qualitatively discuss the performance 
of the lumped-element superconducting oscillators used 
in the measurement, while in section Ild we present a de- 
tailed demonstration of the multiplexed QCM technique. 
In section Ilia, we present a characterization of the two- 
qubit system, where we estimate the qubit parameters 
and coupling energy using microwave spectroscopy and 
an analysis of the ground-state capacitance. In section 
Illb, we use these parameters to estimate the ground- 
state concurrence, and discuss entanglement in the two- 
qubit system. 



II. MULTIPLEXED CAPACITANCE 
MEASUREMENT 

A. Overview 

Let us review the essential details of the quantum ca- 
pacitance measurement for one qubit. For a single paral- 
lel lumped-element LC tank circuit coupled capacitively 
to a single qubit and transmission line, the overall capac- 
itance of the oscillator when the qubit is in its ith energy 
eigenstate is 
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where Ct is the tank circuit capacitance, Cc is the cou- 
pling capacitance between the tank circuit and the 50 Q. 
transmission line, C^j? is the RF gate capacitor, and Cj 
is the combined qubit junction capacitance, Ei is the ith 
eigenvalue of the qubit, and Ug = CgVg/2e is the nor- 
malized qubit gate charge, where Cg is the control gate 
capacitance. The fourth term in Eq. ^ is referred to as 
the quantum capacitance, and is proportional to the cur- 
vature of the qubit energy level. For a single qubit in the 
two- level approximation, where a = Ej/AEc <SC 1, the 
quantum capacitance in the ground and first excited 
(— ) states is given by 
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where C^, is the total capacitance of the qubit island. 
For the coupled two-qubit system, the eigenstates are 
best evaluated numerically, but the quantum capacitance 
will approximately reduce to Eq. ([s]) when one of the two 
qubits is far from its degeneracy point. By measuring the 
phase shift of a reflected RF signal, one can extract the 
quantum capacitance, as described in section lid. Since 
the ground and first excited states have opposite curva- 
ture at the degeneracy point, this technique can be used 
to measure the state of the qubit directly at its operating 
point. 



B. Experimental Design and Setup 

Next we describe a multiplexed version of QCM for 
two charge qubits based on the single Cooper-pair box 
(SCB). Two parallel lumped-element tank circuits with 
different inductances are capacitively coupled to a sin- 
gle transmission line, which is probed with a two-tone 
RF signal. The reflected signal is demodulated in a ho- 
modyne technique with two analog quadrature mixers. 
By measuring the phase shift of each signal, we may si- 
multaneously measure the quantum capacitance of each 
qubit. Such multiplexing schemes have been p revio usly 
employed with the RF single-electron transistor )^^ * ^^! and 
can be readily scaled to read out a large array of qubits 
by using high-Q oscillators. In this experiment, the Q of 
the lumped-element tank circuits is of order 1000. 

The circuit layout and experimental setup is shown in 
Fig. 1. The qubit devices were fabricated using a conven- 
tional shadow-mask evaporation technique. The qubit 
features shown in Fig. lA were patterned with electron- 
beam lithography, and the qubits themselves each consist 
of a pair of small-area Al/AlOx/Al tunnel junctions in 
a DC-SQUID conflguration. Note that the areas of the 
two loops are different, so that the Josephson energies of 
the two qubits can be tuned concurrently with a single 
superconducting magnet. 

The two qubits are coupled by a flxed capacitor, which 
is formed from the shadow of the islands themselves, 
spanned by a small metal bridge. While the bridge ef- 
fectively introduces a pair of tunnel junctions within the 
coupling capacitor itself, the area of these junctions is 
large enough (720 x 230 nm) that the single Cooper-pair 
charging energy of the bridge is negligible. Assuming that 
the junction capacitance scales linearly with the area, we 
estimate this charging energy of the bridge to be 19 mK. 
The tank circuits are coupled directly to the qubit islands 
with large interdigitated capacitors. 

The parallel LC tank circuits, which are shown in 
Fig. IB, are patterned with photolithography and made 
from a superconducting Al/Ti/Au trilayer with respec- 
tive thicknesses of 300, 200, and 200 A. This material 
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FIG. 1: A) SEM image of the qubit structures for sample 
1. Fabrication is performed using electron-beam lithography 
and double-angle evaporation. The interdigitated capacitors 
in the center are the RF gates, while the short leads to the 
sides are the control gates. The qubit islands are the thin 
strips at the center, and include one side of each interdigitated 
capacitor. The two qubits are coupled with a fixed coupling 
capacitor formed by the shadows of the qubit islands, which 
are spanned by a small bridging structure. B) SEM image of 
the optically patterned tank circuits for sample 1. The left 
spiral inductor has 15 turns (Li — 380 nH) while the right has 
17 turns (L2 = 430 nH). C) Circuit diagram showing qubit 
structures and tank circuits. 



acts as a superconductor with a suppressed transition 
temperature Tc — 450 niK,!^ which acts as a trap for 
quasiparticles in the leads. Fig. IC shows the circuit 
layout. Significant high-frequency filtering of the RF sig- 
nal is performed with a CuNi microcoaxial transmission 
line thermally anchored to each refrigeration stage, and 
a cold circulator mounted to the mixing chamber is used 
to isolate the reflected signal return path from the feed- 
line. The resonant frequencies of the right and left tank 
circuits are 545 and 588 MHz, respectively, which are an 
order of magnitude smaller than the qubit energy level 
spacing. 

In Fig. IC, the SCB junctions for qubits 1 and 2 are 
labeled J_ri.2 and Jli,2- The qubit islands are capaci- 
tively coupled to the LC oscillators through the RF gate 
capacitors, labeled Cfipi,2- These RF gates are the large 
interdigitated capacitors shown in Fig. lA. The oscil- 
lators themselves are lumped-element parallel LC tank 
circuits, the components of which are labeled L1.2 and 
Cti,2 in Fig. IC. Finally, the tank circuits are coupled to 
the 50 transmission line through RF coupling capaci- 
tors Cci and Cc2- Nominal tank circuit parameters are 
Li ^ 380 nH, L2 = 430 nH, Cti = 172 fF, Ct2 = 177 
fF, and Cci = Cc2 = 20 fF. 

The quality factor Q of the oscillators are on the order 
of 1000, as opposed to ~ 20 for typical implementations 



of the RF-SET and previous QCM measurements. This 
higher Q value is due to very small coupling capacitances 
between the LC oscillators and the 50 transmission 
line. This results in a larger phase-shift of the reflected 
signal for a given quantum capacitance, thus permitting 
a more sensitive measurement and allowing the use of a 
smaller RF gate capacitor. This reduces the noise cou- 
pled from the RF line onto the qubit, as does the flltering 
effect of a high Q oscillator. Furthermore, a higher Q os- 
cillator permits closer frequency spacing of the two oscil- 
lators, so that both resonances can fit within the 50 MHz 
bandwidth of the cold circulator. When the oscillators 
are monitored with RF reflectometry, the loss of signal 
bandwidth associated with a high Q oscillator limits the 
effectiveness of this circuit for short-pulse or single-shot 
readout. However, when the oscillators are allowed to 
self-oscillate by providing balanced feedback, a change in 
the oscillator frequency can in principle be detected in a 
time much shorter than the oscillator ringdown time.l^^ 
Demodulation is performed by sphtting the amplified RF 
output signal into two quadrature mixers. 



C. Oscillator Performance 

In an ideal superconducting tank circuit, there is min- 
imal intrinsic loss and the quality factor is set simply 
by the coupling capacitance to the 50 f2 transmission 
line. However, in direct measurements of the linear re- 
sponse of the tank circuit as a function of temperature 
and magnetic field, we represent the intrinsic losses of 
the tank circuit at the resonance frequency with a series 
resistance of 1 in devices fabricated on single-crystal 
quartz substrates. This value of the resistance was deter- 
mined by fitting the tank circuit response to a lumped- 
element model. The series resistance manifests itself as a 
decreased amplitude and complex offset in the observed 
signal, which must be carefully accounted for when pro- 
cessing the data. 

For an applied RF power of -130 dBm, the current 
in the inductor is 0.4 nA, and this anomalous resistance 
corresponds to a power dissipation of roughly -160 dBm 
in the inductor. Measurements of the temperature and 
magnetic field dependence of the resonant response of 
similar tank circuits were consistent with predictions of 
a standard two-fluid model for the kinetic inductance and 
surface resistance of the superconducting trilayer. How- 
ever, at low temperature and zero fleld, this model pre- 
dicts negligible dissipation, so this series resistance can- 
not simply be ascribed to the AC surface resistance of the 
superconducting trilayer. Furthermore, an all-aluminum 
tank circuit fabricated with electron-beam lithography 
shows comparable losses. 

Likewise, a straightforward estimate of the power dis- 
sipated into eddy currents in the Au-plated Cu sample 
box yields a dissipated power of -200 dBm, which is not 
enough to explain the loss. Another possible dissipa- 
tion mechanism is electromagnetic loss in the substrate. 



4 



Finite-element simulations predicted a small dissipative 
component in the resonant response given a material loss 
tangent typical of single-crystal quartz, but zero dissipa- 
tion when given a loss tangent typical of sapphire. Sim- 
ilar devices fabricated on polished R-planc sapphire had 
an inductor series resistance of 0.20 ± 0.08 fi. However, 
the high dielectric constant of sapphire required a com- 
plete tank circuit redesign, with significantly less on-chip 
metal, so it is difficult to say with certainty whether the 
improvement is attributable to the improved loss char- 
acteristics of the substrate alone. All data shown in this 
paper was collected from devices using quartz substrates. 
While the tank circuit dissipation mechanism is still un- 
der investigation, significant improvement can be practi- 
cally achieved simply by increasing the coupling capaci- 
tances Cci,2 to the external transmission line. 
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D. Demonstration of Multiplexing 

A demonstration of the multiplexed QCM technique 
using sample 1 is shown in Fig. 2. In this experiment, 
tank circuits 1 and 2 are probed simultaneously, and the 
two qubit control gate voltages are ramped concurrently 
in the same direction. Equal voltage ramps from ±3 mV 
are applied to both qubit control gates, with a ramp fre- 
quency of 104.4 Hz. Note that in this configuration the 
system does not necessarily pass through the two-qubit 
mutual degeneracy point. In this experiment, the sample 
is mounted on the mixing chamber of a dilution refriger- 
ator at its base temeperature of 20 mK. In this analysis, 
the capacitance signal is presumed to be 2e-periodic, i.e. 
there is negligible tunneling of non-equilibrium quasipar- 
ticles. When quasiparticles tunnel across the SCB junc- 
tions, the gate charge switches rapidly by one electron, 
and the time-averaged capacitance signal is a weighted 
average of the shifted and unshifted capacitance peaks. 
Such a signal is typically referred to as "le-periodic" . 
The assumption that the observed signal is 2e-periodic 
at low temperatures is based on the observation of a 
transition to a le-periodic signal between 250 and 300 
mK, which is consistent with the thermal occupation of 
equilibrium quasiparticle states. 

Figs. 2A-B show raw in-phase (I) and quadrature (Q) 
oscilloscope traces for both qubits, as labeled in the figure 
caption. These traces were recorded simultaneously. The 
X-axis displays the sweep time of the qubit gate ramp. 
Figs. 2C-D show the phase shift and magnitude response 
for both qubits, which represent the reactive and dissi- 
pative components of the qubit signal, respectively. As 
discussed above, the qubit response is transformed by the 
loss in the tank circuit, adding an offset in the complex 
plane. To compensate for this effect, the signal is re- 
centered by subtracting this offset before computing the 
magnitude and phase shift. 

For qubit 1, the observed phase shift is 171°, while for 
qubit 2 the phase shift is 117°. The difference in the over- 
all magintude of the phase shifts for the two qubits is due 



FIG. 2: (Color Online) Illustration of multiplexed QCM tech- 
nique. Dark grey (blue) and light grey (red) curves are data 
sets taken simultaneously for qubits 1 and 2, respectively. 

A) Raw oscilloscope traces for in-phase (I) signal compo- 
nent. The X-axis is the sweep time of the qubit gate ramp. 

B) Raw traces for quadrature (Q) components of both qubit 
signals. C) Extracted phase shift for both qubits in degrees, 
accounting for the loss in the tank circuit. D) Signal am- 
plitude for both qubits. E) Quantum capacitance signal for 
qubit 1 extracted from the phase shift, as described in the 
text. The black (smooth) curves are fits to Eq. (|3| with pa- 
rameters listed in the text. F) Quantum capacitance of qubit 
2. 



to the difference in i?ji,2 for this particular value of the 
DC flux bias, which was tuned to maximaize the phase 
shift for qubit 1 . The observed phase shifts are quite large 
compared to previo us QC M experiments, which were typ- 
ically less than 10°. This is a result of the higher Q 
value for the oscillators and the large RF gate capac- 
itances. In figure 2C, the phase response quantifies the 
dispersive response of the oscillator to the changing qubit 
ground state, as discussed in sec. II. The amplitude re- 
sponse, as shown in Fig. 2D, corresponds physically to 
the absorption of RF probe power by the qubit system 
itself. Since the measurement shown in Fig. 2 is per- 
formed with the qubits in their ground state, the gate 
voltage dependence of the signal amplitude is relatively 
weak. For qubit 1, the visibility of the gate voltage de- 
pendence is approximately 7%, while for qubit 2 it is 
below the noise level.--' 

To extract the quantum capacitance from the phase 
shift, we use a lumped-element circuit model with the 
nominal parameters given in Sec. lib. In the case where 
the series resistance due to the intrinsic oscillator loss is 
neglected, the quantum capacitance is given by 
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where uj is the driving frequency, Uc = l/\/LCc, Z = 
(1 + cos (^(rig)) /sin (/>(ng) is the overall impedance 
of the tank circuit, Zq = 50 is the characteristic 
impedance of the transmission line, and 4'{ng) is the 
phase shift of the reflected signal as a function of gate 
voltage. This formula is used to extract the quantum ca- 
pacitance as shown in Fig. 2E. The noise at the peaks is 
due to the fact that Eq. Q is a rapidly varying function 
as 4>{ng) — !■ 0, which occurs at the degeneracy point. 
Note that for clarity, the quantum capacitance traces 
for qubit 2 have been shifted horizontally by 0.5 Cooper 
pairs. The solid curves in Fig. 2E are theoretical plots of 
Eq. (H for Eci/h = Ec2/h = 190 mK, Eji/h = 50 
mK, and Ej2/kb — 120 mK. For this experiment, the flux 
bias was tuned to minimize Eji. These qubit parameters 
are consistent with single-qubit microwave spectroscopy 
measurements, performed when one of the two qubits was 
far from its degeneracy point. The relatively small value 
of Ec is due to the large RF gate capacitors, as shown 
in Fig. lA. 



III. COUPLED QUBITS 

A. Characterization of Coupled Qubits 

We apply the measurement technique discussed above 
to characterize the coupled two-qubit system. We use 
two characterization techniques in two different samples. 
In sample 1, the qubit parameters were estimated by mi- 
crowave spectroscopy. In sample 2, the qubit coupling 
energy was inferred by mapping the ground state capac- 
itance as a function of both qubit gate voltages. Sample 
2 had a similar design to sample 1, but with increased 
charging energies and an increased coupling capacitor. 
Using the qubit parameters for each sample, we estimate 
the concurrence and discuss the ground-state entangle- 
ment. Unfortunately, technical limitations prevented us 
from characterizing both samples using both methods. A 
summary of the extracted parameters is shown in Table 
1. 

The coupled-qubit Hamiltonian is given by Eq. ([T|). 
The coupling between the two qubits introduces avoided 
level crossings in the energy spectrum. The coupling also 
introduces a gate voltage asymmetry into the energy level 
diagram. By probing the system with fixed-frequency 
continuous-wave (cw) microwaves as a function of both 
gate voltages, we can map out the lowest-lying energy 
levels of the two-qubit system. 

Fig. 3 shows data from such an experiment performed 
with sample 1, taken at a fixed microwave frequency of 
11 GHz. The oscillator phase shift of qubit 1 is recorded 
while the system is perturbed by cw microwaves, at a 
range of gate voltage points for the two qubits. The 
phase shift signal taken without microwaves was sub- 
tracted from the data, so the plot in Fig. 3 shows the 
deviation in the phase shift due to the microwave exci- 
tation alone. A voltage ramp is applied to the gate of 




FIG. 3: Mesh plot of the phase shift deviation due to 11 
GHz microwave excitation, as a function of both qubit gate 
voltages, performed using sample 1. Gate charge units are 
in Cooper pairs. Phase shift data taken without microwaves 
has been subtracted to enhance the visibility of the spectro- 
scopic signal. Black lines are numerical calculations of energy 
eigenvalues commensurate with 11 GHz microwaves. Solid 
parallelograms correspond to transitions to the first excited 
state, while the dashed parallelograms correspond to transi- 
tions to the second excited state. The degree to which the 
parallelograms tilt is a good indicator of the coupling energy 
= 25 ± 5 mK. 



qubit 1, while the gate of qubit 2 is held at a fixed po- 
tential. After accumulating a time-averaged phase shift 
signal as demonstrated in Fig. 2, the gate voltage for 
qubit 2 is stepped to the next value. To correct for long- 
time DC drift in the qubit 2 gate charge offset, each av- 
eraged trace is shifted vertically to fix the location of the 
quantum capacitance peaks at half a Cooper pair. This 
also corrects for the effects of cross-capacitance between 
the two qubit gates. Finally, the unperturbed peaks are 
subtracted. The blue represents a large negative phase 
deviation, indicating a more positive quantum capaci- 
tance, while the red is the reverse. The blue parallelo- 
grams correspond to transitions to the first excited state, 
which has a more positive quantum capacitance than the 
ground state. The red spots correspond to transitions 
to the second excited state, which has a more negative 
quantum capacitance. Again, in this analysis the signal 
is assumed to be 2e-periodic because of the observation of 
a transition to le-periodicity between 250 and 300 mK, 
as equilibrium quasiparticle states become occupied. 

The solid lines are numerical plots of the difference 
between the energy levels of Eq. ([T|) , at an energy corre- 
sponding to an 11 GHz microwave excitation. In these 
plots, Eci/h = Ec2/kb = 190 mK, Eji/h = 340 mK, 
Ej2/kb — 430 mK, and Em/kt — 25 mK. These values of 
Ec 1,2 are consistent with those extracted from microwave 
spectroscopy performed on each qubit independently at 
a similar fiux bias. Note that the value of the fiux bias 
is different in this experiment compared to the measure- 



6 



ments discussed in section Ild, leading to larger values 
of -Bji,2- The left-leaning parallelograms indicate tran- 
sitions from the ground to the first excited state, while 
the right-leaning parallelograms indicate transitions from 
the ground to the second excited state. Transitions to 
the third excited state are not expected at a microwave 
frequency of 11 GHz for these values of the qubit pa- 
rameters. The coupling energy Em was extracted largely 
by inspecting the amount of tilt in the parallelograms, 
and the width of the splittings. In the absence of cou- 
pling, the parallelograms revert to rectangles. This tech- 
nique conservatively permits extraction of the coupling 
energy within an error of ±5 mK. For sample 2, which 
had smaller RF gate capacitors and hence a smaller quan- 
tum capacitance, the signal-to-noise ratio was much lower 
and experiments of this type could not be conclusively 
analyzed. 

We also demonstrate a second technique for establish- 
ing the coupling energy. As can be seen from Eq. ([ij, 
the coupling between qubits renormalizes the effective 
charging energies, so that the qubit 1 charging energy 
((Tzi term) depends on 71^2 as well as n^i, and vice versa 
for qubit 2. In a two-dimensional plot of the ground 
state capacitance, this leads to an electrostatic "kink" 
feature in the quantum capacitance at the mutual de- 
generacy point. In effect, an excess charge on the island 
of one qubit changes the potential on the island of the 
other. The depth of this kink depends strongly on the 
interqubit coupling energy, which can be extracted by 
comparing the ground state capacitance as a function of 
both gate voltages to the theoretical values. This experi- 
ment is in some respects similar to previous ground state 
characterizations of three- and four-qubit systems in the 
flux domain.!^ 

The results of such an experiment performed using 
sample 2 are shown in Fig. 4. This sample had much 
smaller RF gate capacitors than used in sample 1, which 
gives a smaller phase shift signal but is less disruptive to 
the qubit state. Decreasing the RF gate capacitance also 
increases the charging energy, which is estimated from 
cw spectroscopy to be Eci/kb = £02/ kt = 740 mK in 
this sample. Sample 2 also had a significantly increased 
coupling capacitance, so that the coupling energy 
was held to the same order of magnitude as in sample 1. 
The data in Fig. 4 is the phase shift in oscillator 1 as a 
function of both gate voltages Ugi and ng2. 

The black line in Fig. 4 is a theoretical calculation of 
the gate charge location of the center of the quantum 
capacitance peak. This calculation was done by numeri- 



Sample 


Eci 


EC2 


Eji 


Ej2 Em 


Method C 


1 


190 


190 


340 


430 25 


Spectroscopy 0.06 


2 


740 


740 


290 


290 80 


Ground State 0.27 



TABLE I: Summary of extracted parameters for samples 1 
and 2. All energy values are scaled by the Boltzmann constant 
and given in mK. The concurrence, C, is unitless. 




Qubit 1 Gate Charge (ngi) 

FIG. 4: (Color Online) Extraction of the coupling energy by 
fitting the ground state capacitance as a function of both 
gate voltages, performed using sample 2. The data is the 
phase shift in degrees for oscillator 1 as a function of both 
gate voltages. The black line is a theoretical plot of the gate 
voltage location of the quantum capacitance peak. From this 
data, we extract a coupling energy of Em = 80 ± 20 mK. 

cally computing the ground state energy of two coupled 
qubits in a sixteen-level approximation (four levels for 
each qubit), and computing the ground state capacitance 
from Eq. In this simulation, the extracted parame- 
ters were Eci/kb = Ec2/h = 740 mK, Eji = Ej2 = 290 
mK, and Em = 80 mK. This technique is somewhat less 
sensitive to the parameter values than the spectroscopic 
technique described above, with an uncertainty in the 
extracted value of the coupling energy of ±20 niK. In 
this plot, the values of Eci,2 and -E,/i,2 are taken from 
cw spectroscopy measurements far from the mutual de- 
generacy point, so the only freely adjustable fitting pa- 
rameter is the coupling energy Em- This technique could 
not be convincingly applied to data taken with sample 1, 
since the qubit parameters were such that the depth of 
the electrostatic excursion was a much smaller fraction 
of the total width of the capacitance peak. As a result, it 
could not be reliably separated from the overall DC drift 
of the gate charge. 



B. Entanglement of Two Qubits 

In studying coupled pairs of qubits in the context 
of quantum information processing, it is important to 
quantify the degree of entanglement between the two 
qubits. A valuable entanglement measure for both pure 
and mixed states of bipartite systems is the concurrence 
C, which is related in a straightforward way to the en- 
tanglement of formation.l^il The concurrence is an entan- 
glement monotone which ranges from for a completely 
separable state to 1 for a maximally entangled state. 

For a pure state {ip), the concurrence of the two-qubit 
system is defined by 
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FIG. 5: (Color Online) Ground-state concurrence of two 
qubits. A) The concurrence, defined by Eq. ([5| as a function 
of both gate voltages for the qubit parameters extracted from 
sample 2, as shown in Fig. 4. Here Eci/kt — Ec2/kt = 740 
mK, Eji/kt = Ej2/kt = 290 mK, and Em/ks = 80 mK. 
B) Concurrence at the mutual degeneracy point in the special 
case of two identical qubits, plotted as a function of coupling 
energy for three different values of Ej. The green (dashed) 
curve shares the same parameter values as the plot in (A). 



C = \ {il;\ay (g) ay\'ijj*)\ 



(5) 



where is the complex conjugate of {ip). Fig. 5A shows 
a numerical calculation of the concurrence as a function 
of both gate voltages when the system is purely in the 
ground state. This calculation was performed using the 
qubit parameters of sample 2, as extracted from the data 
shown in Fig. 4. Note that the concurrence takes on its 
maximum value of 0.27 at the mutual degeneracy point, 
Ugi ~ ng2 = 0.5, and becomes negligible near the edges 
of the plot. Despite the fact that the strength of the cou- 
pling is fixed, the ground state is completely factorizable 
when the system is far from the mutual degeneracy point. 

In the particular case of identical qubits, where Eci = 
Ec'2 and Eji = Ej2, the problem simplifies considerably. 
When expressed in the singlet-triplet basis, the singlet 
and triplet subspaces of Hamiltonian ([T|) are fully decou- 
pled, and the Hamiltonian can be cleanly diagonalized in 
closed form.'^ At the degeneracy point, the ground state 
of the two-qubit system is given by 



/Ej 



This leads to a 



m 



\9) = 



[|00) + |11) + A(|01) + |10))] (6) 



simple expression for the concurrence, C — -^2ip[ 

the zero-temperature limit. Note from this formula that 
C ^ as Em 0, and C ^ 1 as Em oo. This 
expression is plotted in Fig. 5B as a function of coupling 
energy, for three different values of Ej. Note that the 
green (dashed) curve corresponds to the qubit parameters 
used to construct Fig. 5A, Eji/kb — £'72/^6 — 200 mK. 
The blue (solid) and red (dotted) curves correspond to 
Eji/kb = Ej2/kb = 50 mK and Eji/kb = Ej2/kb = 400 
mK, respectively. Using the qubit parameters extracted 
in section II, we find that the concurrence at the mutual 
degeneracy point Ugi = 71^2 = 5 is C = 0.06 for sample 
1 and C = 0.27 for sample 2. Note that at the mutual 
degeneracy point, the concurrence does not depend on 
Eci.2, although it is strongly dependent on Eji^2 and 

Em- 



TV. CONCULSIONS 

We have used a multiplexed quantum capacitance mea- 
surement technique to characterize a system of two en- 
tangled superconducting qubits. We have determined the 
energy scales of the two-qubit system, both through mi- 
crowave spectroscopy and an examination of the ground 
state capacitance. From this information, we have esti- 
mated the concurrence, the degree of entanglement be- 
tween the two qubits, and discussed the ground state 
entanglement for superconducting qubits. 

In the multiplexed QCM technique, two on-chip su- 
perconducting lumped-element LC oscillators are capac- 
itively coupled to the qubit islands, and monitored thor- 
ough a single RF line. Since the overall capacitance is 
dependent on the qubit state, a measurement of the os- 
cillator with RF reflectometry constitutes a dispersive 
measurement of the qubit, which is applicable directly 
at the degeneracy point. This technique is readily scal- 
able to read out a large array of qubits, which is a key 
ingredient in the development of a large-scale quantum 
computer. 
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